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Abstract 

We study spherically-symmetric structures in Conformal Gravity and in a scalar-tensor 
extension and gain some more insight about these gravitational theories. In both cases we 
analyze solutions in two systems: perfect fluid solutions and boson stars of a self-interacting 
complex scalar field. In the purely tensorial (original) theory we find in a certain domain 
of parameter space finite mass solutions with a linear gravitational potential but without a 
Newtonian contribution. The scalar-tensor theory exhibits a very rich structure of solutions 
whose main properties are discussed. Among them, solutions with a finite radial extension, 
open solutions with a linear potential and logarithmic modifications and also a (scalar- 
tensor) gravitational soliton. This may also be viewed as a static self-gravitating boson star 
in purely tensorial Conformal Gravity. 

1 Introduction 

Conformal Gravity p] (CG) was proposed as a possible alternative to Einstein gravity ("GR"), 
which may supply the proper framework for a solution to some of the most annoying problems 
of theoretical physics like those of the cosmological constant, the dark matter and the dark 
energy. 

It is therefore very much required to investigate its predictions and consequences as further 
as possible. Here we choose to concentrate in localized solutions and to start an investigation 
of their properties. We take two simple matter sources: perfect fluid and complex scalar field, 
we find localized solutions for both kinds of sources and present their main features. 

The main ingredient of CG is the replacement of the Einstein-Hilbert action with the Weyl 
action based on the Weyl (or conformal) tensor C K \^ v defined as the totally traceless part of 
the Riemann tensor (we use R K \^ V = d v F% — d^F^ + ...): 

1 R 

CkX^u = Rk\/j,v - -^(gKfiRxu - 9kvR\h + 9\vRkh - gXpiRnv) + —{9kij,9Xu ~ 9kv9\p), (1.1) 



so the gravitational Lagrangian is 



£g — ~ ~2^f^ K ^ v ^ ^ (1-2) 
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where a is a dimensionless positive parameter. The gravitational field equations are formally 
similar to Einstein equations where the source is the energy-momentum tensor and in the 
left-hand-side Bach tensor W^ u replaces the Einstein tensor: 



W, 



—T 
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(1.3) 



Bach tensor is defined by: 



W, 



1 



V^V.i? - V A V A iV + ^(R 2 + V A V A i? - 3R kX R kX ) 9 ^ + 2R KA R^ uX - =-RR^ (1.4) 
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Since Bach tensor is traceless, the energy-momentum tensor must "comply" so we will consider 
only sources with = 0. 



The general spherically-symmetric line-element may be simplified by exploiting the confor- 
mal symmetry and has the form [I]: 



ds z 



B(r)dt 2 - dr 2 /B(r) - r 2 (d6 2 + sin 2 



The non- vanishing components of Ricci tensor and the Ricci scalar are 
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A useful property of these components is the following: 

B(rB)"" 



W% - W r r 
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In vacuum this is easily integrated to give 

B(r) = cq + C\r + C2/V + nr 2 
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(1.11) 



where the relation between the coefficients comes from the WjT = equation which is of a third 
order. In a non-relativistic fourth order gravity a similar situation is encountered, namely the 
fourth order "Poisson equation" 

V 2 V 2 u = -/i (1.12) 
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where h(r) is the source term. In the spherically symmetric case V 2 V 2 tt = (ru)""/r and u(r) 
is given also by (jl.lip without any relation between the parameters. On the other hand, the 
parameters are related to the source (assumed to extend within r < a) by 

1 f a 1 f a 

c\ = — I r 2 h(r)dr ; C2 = — / r 4 h(r)dr. (1-13) 
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Since k is not fixed by the source, the Kr 2 term may be considered as a possible background 
field or in the relativistic context, a cosmological constant contribution^]. Note also that the 
volume integral of the matter density (i.e. of h{r)) turns up as the coefficient of the linear 
term in the potential rather than the 1/r one. It is related to the fact that in this theory the 
potential of a point particle is linear in accord with the behavior of the Green function. This 
linear potential enables one to explain galaxy rotation curves without assuming dark matter 

mm- 

For the general case of extended sources, we note that since the field equation is of fourth 
order, a special care should be taken with the boundary conditions. It is easy to see that 
u'(0) and u"'(0) should vanish. The value of u"(0) or u"(oo) may be free if solutions with a 
"cosmological" Kr 2 term are allowed. If on the other hand the background is assumed to be 
empty ("flat"), we may impose further u"(oo) = as well. If the source is localized, the second 
derivative at the origin is related to the first moment of the matter distribution as 

1 r a 

u"(0) = - rh(r)dr (1.14) 
3 Jo 

Now let us return to the relativistic field equations with a perfect fluid source described by 
T]f = diag(p,—P r ,—P±,—P±) (with the additional conformal condition = 0). Thanks to 
()1.10p they reduce to a single very simple field equation: 

= <!•»> 

which has a similar structure to the fourth order Poisson equation (|1.12p . By comparison we 
notice that taking a > corresponds to gravitational attraction in the weak field limit. 
Equation (|1.15|) should be solved together with the conservation equation 

p; + -(3P r -p) + ^-(p + P r ) = (1.16) 
r IB 

and an additional equation of state which relates algebraically p, P r and P± . Regularity of the 
Bach tensor at the origin introduces an additional boundary condition, -B(O) = 1 to those of 
the Poisson case: B'(0) = B"'(fi) = 0, B"(oo) = 2k. 

The inertial mass of such a spherical solution is the ordinary 

/poo 
d?xy/\g\ T ° = 4vr y r 2 p(r)dr. (1.17) 

However, since the potential of a point particle in this theory is linear, the gravitational mass 
is identified as the coefficient of the linear term in the vacuum potential - see (|1.13|) and (|1.15|) : 



Mq = 12tt / drr 2 (p(r) + PJr))/B(r) = a 

Jo a 



(1.18) 



to be concrete, R — 4A = —12k, so k > corresponds to AdS. 
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The other parameter, C2 (the coefficient of the 1/r term in the potential) has a dimension of 
length which by utilizing Newton's constant can be converted to a mass. However, we don't 
have an appropriate dimensionful parameter at our disposal, so we will call C2 the "second mass 
parameter" . In terms of the source functions it is given by 

c 2 = - drr 4 (p(r) + P r {r))/B(r) (1.19) 
4 Jo 

We will see in the next sections that this integral is not always convergent, and whenever it 
does, it has the wrong sign for an attractive force, causing a non-Newtonian "near field" of such 
sources. Actually, this problem that ordinary continuous sources do not produce a Newtonian 
component in CG was noted already by Mannheim and Kazanas [3] (following even earlier 
studies [H El [7] from the 1960's and 70's). Mannheim and Kazanas pointed out towards a 
possible solution based on the fact that a highly singular source can produce a potential with 
both ci > and c% < 0. Still, when the implications and consequences of CG are analyzed, 
smooth matter distributions should be considered and studied since they are more widely used 
to model astrophysical and cosmological sources. 

2 Spherically- Symmetric Perfect Fluid Solutions 

In accord with our objective, which is investigating the properties of self gravitating solutions 
in CG, we solved Eqs. (|1.15p . (|1.16p for a set of matter distributions. 

The simplest of all sources is a constant energy density, p(r) = po (for r < a and outside), 
but unlike the Einsteinian case, there are no finite mass solutions of this kind in our case. 

The "next to simplest" source is a polytrope - either linear with P r = p/noi non-linear (and 
anisotropic) with P r = Pq(p/2>Pq)~ i where n, 7 and Po are all positive constants. The parameter 
Po is indeed the central value of the pressure (if P r (0) is finite). Note that the special value 
7 = 1 gives only the n = 3 case of the linear relation which corresponds to isotropic radiation. 
The other values of n cannot be obtained as a limit of the non-linear polytrope. 

Next we move to general polytropes, that is, density and pressure related by 

p = 3P (P r /P ) 1+A (2.1) 

where for convenience we parametrize the polytropic index by I/7 = 1 + A. The construction 
of regular solutions for r £ [0, 00] requires the boundary conditions 

P(0) = 1 , P'(0) = , B"'(0) = , P r (0) = Po , (2.2) 

the fifth boundary condition was fixed by imposing the value P"(oo) which is related to the 
free "cosmological constant parameter" k (see (jl.lip ). The numerical results further indicate 
that the solutions behave asymptotically according to 

B = Kr 2 + B x r + P + • • • , P- oc (2.3) 

where the constant p depends on k and on A. 

We will discuss separately the solutions available for vanishing and non vanishing n, that 
is P"(oo) = and P"(oo) ^ 0. 
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Solutions with k = 



By examining the conservation equation (|1.16p . we obtain the physically acceptable decay 
of the function P r (r) in terms of the parameter A. It turns out that solution with an asymptot- 
ically decreasing P r can only occur for A > 0. We then get P r ~ r~ 2 for A = and P r ~ r~ 7 / 2 
for A > 0. From these observations, it follows that the inertial and gravitating mass given 
by (|1.17|) and (|1.18|) . do not converge in the case ^4 = ("radiation ball"). Constant density 
solutions (A = — 1) do not exist as well. 

Expanding Eq. (j!.16[) around the origin, we further observe the following relation 

3 P"(C\) 

( l_^)£m + 2 S"(0) = (2.4) 

suggesting that the value A = 2/3 should play a role in the solutions. Integrating the equations 
we obtained finite mass solutions for < A < 2/3. Fig. contains graphic representations of 
three solutions in this range. Actually, we solved a dimensionless version of Eqs. (|1.15p . (|1.16p 
and (|2.ip for B, P r /Po and p/Po in terms of x = r(aPo) 1//4 - It is clear from the plots that the 
gravitational potential is asymptotically linear which is the required form in order to explain 
the galactic rotation curves within this context j2[ [T]. However, a closer inspection shows 
that the 1/r component, which is necessary for the recovery of the Newtonian (Schwarzschild) 
behavior in smaller scales, is missing. This is reflected by the fact that the coefficient of the 
1/r term, C2 (see (|1.19p ) diverges. 

Several parameters characterizing the solutions (namely the masses, the values B"(0), P"(0) 
and .B'(oo) are depicted on Fig. [lb; it shows in particular that the solution is well defined in 
the limit A = 0. In fact, in this case, we have P r = Pq/B 2 but the masses are infinite. The 
limit A — > 2/3 is more subtle. It seems indeed that in this limit the function B(r) approaches 
B = 1 on the full space, while the function P r becomes more and more concentrated around 
the origin and |P"(0)| — > oo. We checked that the relation (|2.4p is obeyed. At the same time 
the masses approach zero. 

It is expected to also consider the equations for A > 2/3. We were able to obtain solutions 
in this case. Our numerical results strongly suggest, however that no globally regular solutions 
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(a) 



(b) 



Figure 2: Perfect fluid solutions with n = 0.1: (a) the profiles of three solutions; (b) plots of several 
characteristics of the solutions as a function of the parameter A. Note that unlike the k — case, here 
B"(0) does not vanish as A — ► 2/3. 

exist there. The solutions have B"(0) < and the function B(r) approaches a zero at some 
finite value r = tq. At the same time, the pressure becomes singular for r — > ro, suggesting 
that the solution is singular. Of course, the numerical construction of such solutions cannot 
be achieved directly for r G [0, oo]; in fact, we proceed on a small interval r £ [0, r max ] and 
gradually increase r max . 

Solutions with k > 

The pattern of the solutions is very similar to the case k = 0. In particular, regular solutions 
are also limited to A < 2/3 and singularities appear for A > 2/3. It is worth noticing that the 
conservation equation (|1.16j) implies now P r ~ r~ 4 (for A > 0), so that the masses are finite in 
the limit A — > 0. Of course the values of the masses depend on the value adopted for B"(0). 
Another difference with respect to the case k = 0, is that now the second mass parameter 
(Eq. (|1.19p ) converges and a 1/r term appears in the gravitational potential, but with a wrong 
sign. Moreover, the convergence is related to the non-vanishing cosmological term k, so this 
potential is of a universal nature rather than of local one. 

In the case k < 0, the field B(r) has a node at a finite r say, r = tq which leads to a 
singularity of the matter function. This is just the de Sitter horizon which is related to the fact 
that de Sitter space does not admit a globally static coordinate system. We will not consider 
the possibility of asymptotically de Sitter space further in this work apart from few mentions. 

3 Boson Stars 

Among all the higher order gravitational theories [HI [9], CG is unique in the sense that it is 
based on an additional symmetry principle. The conformal symmetry imposes severe limitations 
on the allowed matter sources. When matter is described in terms of an energy-momentum 
tensor it should be traceless as mentioned above already. Similarly the matter Lagrangian 
is very much constrained, but the Abelian Higgs model is essentially still consistent with the 
conformal symmetry provided the scalar field "mass term" is replaced with the appropriate 



"conformal coupling" term which introduces a coupling to the Ricci scalar R. The matter 
Lagrangian which we will use here is therefore 

C m = l -(D^T(D^) - li?|$| 2 -\\<S>\ 4 - \F, V F»\ (3.1) 
and the resulting field equations are 

D ll D l *$ + \\<&\ 2 $ + —$ = (3.2) 

IP 

VpF*" = ~[$*(D V $)-$(D V $)*] = J v . (3.3) 
The gravitational field equations are (|1 .3j) with 

V = T^ inima + l - (<^V A V A |<1>| 2 - V M V,|<1>| 2 ~ G^\ 2 ) (3.4) 

Tftu being the ordinary ( "minimal" ) energy-momentum tensor of the Abelian Higgs model 

and G^y is the Einstein tensor. 

The simplest spherically-symmetric localized solution of this system is the boson star [10} 
HU [12] which requires a global U(l) symmetry only - that is = and <3? = f(r)e lult . This 
yields a global conserved charge which is responsible for its existence. 

The components of the energy-momentum tensor are (after use of the ^-equation (|3.2p ): 

T o = ^^ll + E f 2_ A f + E { fy + 1( B " + ^ + f (3.5) 

6 B 6 12 J 12 w ; 18 1 ■■ •■- 1 
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Since there is only one independent metric component, it is obvious that not all the field 
equations (jl.3p are independent. Actually there is only one independent equation and we may 
use the third order one 



w: 



-T T 



0. 



(3.8) 



However, a much simpler form is again obtained by using (jl.lOp giving therefore the following 
fourth order equation for the metric component B(r): 
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For the scalar field we have the second order equation 
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(3.9) 



(3.10) 



where one should also write explicitly R = 2(1 — B)/r 2 — AB' jr — B" by (|1.6p . 

The inertial mass and gravitational mass of these boson stars are given by equations like 
(|1.17p and (j!.18j) with the necessary adaptations: 



Mi = 4vr / drr 2 T$(r) 



(3.11) 
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;>oo 

M G = l2ir drr 2 {T^{r)-T^{r))/B(r) (3.12) 

JO 

where Tg and T£ are given by (|3.5p and (|3.6p . The second mass parameter is defined in analogy 
with ([LT9jk 

c 2 = f ^ drr 4 (T °(r) - ^(r))/B(r) (3.13) 
The boson star has also a global charge (particle number) which is given by 

poo 

Q = 4ttuj / drr 2 f 2 {r)/B(r). (3.14) 
J o 

It is interesting to note that the field equations form an autonomous system as a result of 
the transformation 

B(r)=V(r)r 2 ; f{r)=ip{r)/r ; r = 1/u (3.15) 
and they also simplify considerably: 

V W + (^) 2 -^" + ^ = (3.16) 

(vy)' + ^+^^-A^ = o. 

Here of course ' = d/du. These field equations may be obtained from the following "reduced 
Lagrangian" : 

There is also a "conserved energy" K (such that K' = 0) : 

* = I «0 2 - ZV'V'") + |(^') 2 + ^ + - ^ 2 - ~y (3.18) 

whose value is not free but fixed to be K = 2/3ct since Eq. (|3. 18j) is equivalent to (|3.8p . 
Moreover, if we define a third degree of freedom W, such that W = V", the equations of 
motion can be derived from the following "ordinary" second order Lagrangian 

L 2 = \^'f + + \^ + \VW ~^(W 2 + 2V'W) (3.19) 

4 Boson Stars: Numerical Results 

In absence of explicit solutions (not even to the simple autonomous system), we approached 
the system of equations (I3.9p . (|3.10p numerically. Using an appropriate rescaling r — > Cr and 
/ — ► Ff, the coupling constants A, a scale with a factor C 2 F 2 while u scales by C. Using these 
rescaling, we can set u = a = 1 in the equations and study the solutions for several values of 
the coupling constant A. If we denote by f{x) and B{x) the solution with u = a = 1 and a 
given A, the solutions with general values of uj and a and self-coupling a A are 

/(r) = -£=/M ; B(r) = B(ur). (4.1) 
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It is also easy to see that the charge Q is independent of the parameter oj and the mass scales 
like u/a. 

Since we chose to solve the fourth order equation (|3.9p . Eq. (|3.8p serves as a constraint. 
Taking the derivatives of the left-hand side of (|3,8p with respect to r and eliminating the 
maximal derivatives B"" and /" by using (|3.9p - (|3.10p . leads to an expression which vanishes 
identically. This implies that the two equations we chose to solve guarantee that the combi- 
nation WjT — § T,T is constant so the constraint will be automatically fulfilled (for any r) by a 
consistent choice of the boundary conditions (such that the constant is 0). 

We first discuss the solutions in the case where the function B{r) is asymptotically linear; 
that is to say B"(oo) = 0, or k = in (jl.lip . The relevant set of boundary conditions for 
solutions of this type is 

5(0) = 1 , B'(0) = , B'"(0) = , /'(0) = , £"(oc) = , /(oo) = (4.2) 

For a better understanding of the numerical results, it is instructive to analyze the asymptotic 
possible behavior of the solutions. The asymptotic form of the B field, i.e. B(r) ~ B\r, enforces 
the function f(r) to obey asymptotically an hyper geometric equation whose solutions are of 
the form 

F . /log(wr) \ 
/( r ) = — sm (— g" — + ^1 , r^oo (4.3) 

where Fo,ip are constants. As a consequence, the function f(r) oscillates asymptotically and 
necessarily develops nodes, rendering the numerical integration technically difficult. We manage 
however to construct the solution by replacing the condition /(oo) = by /(ro) = imposing 
by hand the first zero ro of the function f(r). Proceeding this way, we obtained strong numerical 
evidences that a continuum family of solutions exist, labelled by ro. In particular, the values 
B"(0), /(0), Bi are fixed by ro. Unfortunately, the integrated energy density and particle 
number densities of these solutions behave according to 

J dr- (sin (log(ur)/Bi + if)) 2 ~ J dy sin 2 (y + if) ,y = \og(ujr) . (4.4) 

The corresponding mass and particle number are then infinite. In other words, there do not 
exist in this theory boson star solutions with a linear gravitational potential. 

On the other hand, boson stars exist with a quadratic gravitational potential which corre- 
sponds to an asymptotically anti de Sitter space (a negative cosmological constant). Setting 
B"(oo) = 2k, we obtain the asymptotic form 

B(r) = Kr 2 + B ir + B + ... , / = /i/r + / 2 /r 2 + . . . (4.5) 

Finite mass solution needs to impose the stronger decay at infinity such that f\ = 0. We 
obtained strong numerical evidences that such solutions exist, that is we solved numerically the 
field equations for a wide range of the self-coupling parameter A. A typical profile is presented 
in Fig. [3h. with A = 1 and k = 0.1. Several physical characteristics of the solutions are 
plotted for A G [0, 2] in Fig. [3b- Surprisingly, the solutions seems to persist in the absence of 
self-interaction (A = 0). 

Due to the asymptotic behavior mentioned above, the second mass parameter (I3.13P is also 
finite, namely the gravitational potential contains in this case too a wrong sign 1/r term. 
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Figure 3: Boson stars with n = 0.1. (a) the profiles of a typical conformal boson star solution with 
A = 1; (b) plots of several characteristics of the solutions as a function of A: The inertial mass Mj and 
particle number Q are plotted in units 8n while the gravitational mass Mq is given in units Air. 



5 Scalar- Tensor Conformal Gravity 

CG has been criticized from several aspects both phenomenological and formal. Several authors 
claim that predictions in the weak field limit disagree with solar system observations [13] , yield 
wrong light deflection [H] (see however suggestions |15[ [16] for circumventing the difficulties) 
or more generally, the exterior solution (jl.lip with k = 0, c\ > and C2 < (which yields the 
desired behavior) cannot be matched to any source with a "reasonable" mass distribution [17J. 
In this respect we have found in the previous sections that boson stars indeed cannot produce 
such a behavior. On the other hand, the "anisotropic" polytropes (Eq. (12. ip ) may present a 
linear potential, but the 1/r component is missing. 

Other authors find evidence for tachyons or ghosts [18J or raise the fact that only null 
geodesies are physically meaningful in this theory since the "standard" point particle La- 
grangian is not conformally- invariant |19j . 

This last point (and possibly some of the former) can be easily corrected and can serve as 
a starting point for a consistent conformal theory by adding a real scalar field and turning the 
theory into a scalar-tensor theory. The conformally-invariant point particle Lagrangian will be 

L pp = -Sy/g^ v x^x u (5.1) 

where S is a real scalar field with the usual conformal transformation laws. The gravitational 
Lagrangian (|1.2p will be modified to 

Cg = \ {-\ C ^ CKX>MV + \ V ^ S ~ J2- RS2 ~ 1 S4 ) ( 5 - 2 ) 

where v is a possible self-coupling parameter. 

The field equations will be modified accordingly. First of all, there will be an additional 
scalar field equation: 

V^S + i/S 3 + ^5 = 0. (5.3) 
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Second, we turn to the tensorial equations (|1.3p . Technically, the modification is just an 
additional energy-momentum tensor S^ u in the right hand side of (jl.3p namely 

W, u = + l -S^ (5.4) 



where 



(5.5) 



But in principle the scalar field should be considered as a gravitational degree of freedom which 
is stressed by the absence of the coupling constant a in front of S^. 

The simplest case to be studied is static spherically-symmetric vacuum solutions which 
within this framework are obtained by solving the following simplified version of Eqs. (|3.9|i - 
(IBTTUj) for B{r) and S{r): 



(rB)"" 1 
r + B 



BS' 2 - -S 4 + - ( B' + — ) (S 2 )' - -S 2 
2 4 V r y ' 12 







(r 2 



—EpL - - uS s = (5.6) 

where as before, R should be expressed in terms of B(r) using (jl.6p . Since the equation 

W r r = (l/2)5 r r (5.7) 

is of third order, it will be necessary to assure its validity and this will be done as before by 
using consistent boundary conditions. 

The difference with respect to the boson stars discussed above, is that now we may allow 
singular solutions in analogy with the Schwarzschild solution of standard GR. The no-hair 
theorem which precludes black holes with scalar hair is evidently not applicable in the present 
context. 

Actually, one may prefer to study the system in a different gauge where by conformal 
transformation the scalar field is a constant, S(x tM ) = So- This simplifies considerably the 
general field equations (|5.3p - (|5.5p and gives immediately the result R = —6vS 2 . However, after 
transforming to a constant S(x^), one cannot use the "Mannheim gauge" (Eq. (II. 5j) ) any more. 
The metric tensor will have two independent components and the relatively simple expressions 
for Bach tensor Wv will become quite cumbersome. 

We therefore chose to stick to the "Mannheim gauge" and to use S as a second degree 
of freedom. On the other hand the "effective metric" that a point particle experiences is 
9fiv = S 2 g^ u - see (|5.ip . Consequently, the interpretation of the solutions is now quite different: 
it is now g^y which has the physical significance, and the question of the gravitational potential 
should be answered by analyzing <7oo = S 2 (r)B(r) rather than B(r). 

As for the purely tensorial case with a scalar field, the vacuum scalar-tensor theory yields an 
autonomous system as well. We repeat the transformation (|3.15p now with S(u) = S(l/u)/u 
and get the equations of motion 

V"" + (S') 2 - ^ss" = 

(VEV + - uT, 3 = (5.8) 

6 
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the "reduced Lagrangian" 



Lst = \{V"f + ^(E') 2 + ~ (l - ^) E 2 + ^E 4 <r,<„ 



and the "conserved energy" K$t (now -Kst = 2/3): 

Kst = I {(V"f - 2V'V"') + ^(E') 2 + y EE' - is 2 - ^S 4 . (5.10) 
The second order Lagrangian is in this case just (|3.19p without the uj term: 

W = ^(E') 2 + ^E 2 + ^S 4 + VeE' - \{W 2 + 2V'W) (5.11) 
2 o 4 o b 

6 Vacuum Solutions 

6.1 Schwarzschild-Like Solutions 

Equations (|5.6p possess a three-parameter family of explicit solutions given by 



B { r ) = {i+ r/af-^^^ ]+ ^^-^j^%) , S(r) = * (6.!) 

r (1+r/j/a) 2 \r£ r (l+r^/a)^/ 1 + r/a 

where r^, a, Sq are free parameters; o and So fix respectively the scale of the radial coordinate 
r and of the scalar function S(r). This solution describes a family of black hole space-times 
with a regular horizon at r = satisfying B(rh) = 0. 

Assuming for simplicity v = 0, these solutions have the following behavior: 

B'(r h ) = - + - , B(r^0) = — [ h , B(r -» oo) = ^ +0(r) (6.2) 

The metric g^,, is therefore asymptotically anti-de Sitter (or de Sitter) space with a cosmological 
constant A = — 3/(a(a + r^)). These formulas can be generalized for v 7^ but they become 
more involved; the case v = is sufficient to illustrate our results. Fixing the parameters 
and a (or A) , the solutions (|6.ip can be of three different forms according to the value of a : 

• a > : a single horizon at r = r^. 

• — rh < a < : & regular horizon at r = hidden by a doubly degenerate horizon at 
r = f/j > rh with fft = + \J r f l + 16ko 4 . A < in this case. 

• a < — : a regular horizon at r = r/, and, inside, a doubly degenerate horizon at 



r 



^/i < r h with Vh = r h~ \J r h~^~ 16ko 4 . A > in this case. 



Note however that although -B(r) is quite similar to the solutions of the purely tensorial CG, 
the actual gravitational field "felt" by a point particle is very different, since the components 
of the relevant metric tensor are those of = S 2 g fll/ namely S 2 (r)B(r), S 2 (r)/B(r) and 
r S (r). Already here we can notice that S 2 (r)B(r) increases with r much less steeply, and 
actually goes asymptotically to a constant. Moreover, the circumferential radius rS(r) is also 
bounded. These solutions are therefore closed. On the other hand the limit a — > 00 gives rise 
to yet another kind of solution with constant S(r) and a "purely Schwarzschild" B(r). It is 
just a special case of a whole family of open solutions that will be discussed below in sec. 16.2.21 
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6.2 General Black Hole Solutions 



The family of solutions discussed in the previous section are entirely determined by the scale of 
the scalar field and by the value of the horizon and the parameter o. In particular, the values 
of the horizon, of the derivative B'(rh) and of the cosmological constant are not independent. 
However, since the equation determining the metric field is of the fourth-order, more general 
solutions are expected. In absence of a generalization of the explicit form (|6. 1 1) . we investigated 
the equations by numerical methods. The first step in this direction consists of establishing 
the most general set of appropriate boundary conditions. Prior to this step, the following scale 
invariance of Eqs. (|5.6j) has to be fixed : 

r^Cr , S^^ , B^B (6.3) 

where C is a constant. We will fix this arbitrary scale by imposing a particular value for 
Sh = S(rh). So we define a dimensionless scalar fielcH S/Sh and a radial variable x = r/\a\. 
For the vacuum solution (16, ip . this scale fixing yields the relation So = 5^(1 + Xh)- 
Solutions presenting a regular horizon at x = Xh require the following conditions : 

B(x h ) = , B'[x h ) = b , G\ x=Xh = , H\ x=Xh = , S[x h ) = 1 , 5"(oo) = B 2 = 2na 2 (6.4) 

where the symbols Q, 7i represent (respectively) the conditions of regularity of Eq. (|5.6p at the 
horizon and the constraint (|5.T[) : 

Q = 6B'S' - S{\ - — - B") - 6zaB 3 (6.5) 

x h Xh 

H = 2x 2 h {A{B') 2 - S 2 ) + 2x 3 h (B'S 2 - AB'B") + x 4 h {2{B") 2 - AB'B'" + 2B'SS' - 3vS*) - 8 (6.6) 

The normalization chosen for the field S(x) in (|6.4p fixes the rescaling f|6.3[) . The constants 
b, B2 are a priori independent. They encode the deviation with respect to the vacuum solution 
(|6.ip where the relation between them is fixed to give 



2(bx h - 1 



,2 



B 2 = ^3-^ , (6.7) 



as found by eliminating the parameter a from B'(rh) and B"(oo) of (|6.2|) . This demonstrates 
in particular that, to any positive value of B2 (i.e. negative A), two solutions of the form 
(|6.ip are available. One of these solutions presents a doubly degenerate horizon at x^ = 



Xh + \Joc 2 +8/B 2 , corresponding to a < in (|6.ip . 

Our numerical results show strong evidence that the analytic solutions can be deformed for 
generic values of b, B 2 or, put differently, are just special cases of a much wider family of vacuum 
solutions of the scalar-tensor conformal theory. These new solutions can be characterized by 
their expansion around the horizon , 

B{x) = b{x-x h ) + ^-(x-x h ) 2 + ^-(x-x h f + ... , S(x) = l + s(x-x h ) + 0((x-x h ) 2 ) (6.8) 
2 6 

as well as by their asymptotic behavior 

B(x) = — x 2 + B lX + B + — + 0{x- 2 ) , S= — + 0(x- 2 ) (6.9) 
2 x x 



We will still use S for it too. 
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where the parameters B2, 6 have to be set in the boundary conditions while the parameters 
Bi,Bq, B_ and S\ of the scalar field can be determined from the numerical solutions. 

Due to the decay S ~ 1/x, the combination BS 2 (x) approaches asymptotically the constant 
B2S1/2. This is encouraging since it yields a non-degenerate point particle Lagrangian (see 
Eq. (|5,ip ) in the asymptotic region. 

For B2 = 0, the asymptotic expansion involves 'log'-terms, in particular S(x) ~ (5i + 
S2 log(x))/x + 0(l/x 2 ) and the expansion of B(x) is more involved. Black hole solution ap- 
proaching a de Sitter space-time asymptotically (i.e. with negative B2 or positive A) can also 
be constructed, presenting a cosmological horizon at some radius x = x c with x c > x^. 

We now discuss the new solutions for B2 > 0. 

6.2.1 Case v = 

We first discuss solutions in the case v = 0. Two such solutions are presented in Fig. [4] for 
Xh = 0.5, B2 = 2/(1 +x/J = 4/3. Here the analytic solution corresponding to b = 3 is compared 
to the numerical solution corresponding to b = 1. 

A natural question consists of determining the domain of existence of the solutions with a 
fixed Xh in the plane b, B2- Our numerical investigations reveal that, for fixed x^ and B2, black 
holes exist for < b < b max where the maximal value b max depend on Xh,B2- In the limit 
b — > 0, the horizon becomes extremal. For Xh = 0.5, we find b max « 3.7,3.06,2.3 respectively 
for B2 = 7/3, 4/3, 0. The values of b corresponding to the analytic solutions are 

b « {1.17 , 3.41} for B 2 = 7/3 , b = {- , 3} for B 2 = 4/3 , 6 = 2 for B 2 = (6.10) 

3 

The numerical solutions therefore exist for larger values of the parameter 6 than the analytic 
ones. The following table summarizes these results. 





B 2 = 7/3 


B 2 = 4/3 


B 2 = 


bmax 


3.7 


3.06 


2.3 


banalyticl 


3.41 


3 


2 


banalyticl 


1.17 


4/3 


2 



Table 6.1: Summary of results for the parameter b for three values of B^. 

The parameters 62,63,-Bi are plotted as functions of 6 in Fig. [5] for B2 = 4/3 (branches 
labelled '1'). The evolution of the parameters 62, 63 clearly determines the critical phenomenon 
stopping the solution at 6 = b max . The property that solutions do not exist for b > b max 
suggests that a new branch of solutions should occur for 6 < b max , joining the first branch in 
the limit 6 — > b max . This was confirmed by the numerics: we indeed managed to construct 
a second family of solutions presenting this property. The corresponding data is presented 
in Fig. [5] by the lines labelled with a symbol '2'. Decreasing the parameter 6 along the 
second branch, we observe very peculiar properties. In particular the functions S(x) and B(r) 
stop to be monotonically decreasing, but present respectively a local minimum and a local 
maximum at two different radii which are rather close to the horizon. For 6^0, the position 
of the local extrema slowly move to the horizon and result in large variation of the derivatives 
S'(x), B"(x), B"'(x) in the region of the horizon. The numerical results suggest strongly that 
the solutions tend to a configuration where S(x) presents a singularity at the horizon. This 
appears on Fig. [5] where the parameters 62, 63, s are plotted as functions of 6. 

Profiles of three solutions of the second branch are shown in Fig. [6l Note that the effective 
metric exhibits a similar behavior than above and very different from purely tensorial CG: the 
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12 3 4 

X 

Figure 4: First branch vacuum solutions for the case v = 0: Two profiles for B2 = 4/3 with b = 1 and 
b = 3 (analytic solution). Color online distinguishes between the curves. In a B&W version notice that 
5(x)-curves are the two decreasing ones, while those of B{x)S 2 {x) start on the ir-axis. 

gravitational potential which is encoded in S 2 (r)B(r) increases much less steeply and tends 
asymptotically to a constant, the spacetime seems to have only a bounded extension since the 
circumferential radial distance rS(r) has a finite limit as r — » 00 as well as the proper radial 
distance J drS(r)/ y/B(r). 

6.2.2 Case v ± 

For v 7^ 0, the analytic solutions with fixed Xh,B 2 are real as long as the condition 

B\x\ + 6B 2 vx 2 h + 8B 2 - 3u 2 x 2 h -8v>0 (6.11) 

holds. This defines bounds of the parameter v. We have also tried to deform the numerical 
solutions available for v = to the case v 7^ 0. The features of the solutions are basically 
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Figure 5: The two branch structure of the v — solutions. The bullets indicate the corresponding 
analytic solutions. For the meaning of the various parameters, see text. 
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(a) (b) 

Figure 6: Second branch vacuum solutions for the case v = 0: Three profiles for B^ = 1/3 with 6 = 2, 
6=1 and 6 = 0.25. (a) curves for B(x) and B'(x); (b) curves for S(x), xS(x) and S 2 (x)B(x). 



similar. Keeping the parameters b, B2 fixed and increasing u, it turns out that the coefficient 
Si of the scalar field (defined in Eq. (|6.9p ) increases rapidly and diverges when the coupling 
constant v approaches a critical value. For example, setting b = 1, B2 = 4/3, we find that the 
main branch exists for < v < 1.65 while the second branch exists for < v < 0.075. 

Contrary to our expectation, the two solutions available for v = do not converge to a 
common solution while increasing v gradually. 

All the above mentioned solutions are closed ones. However, there exists another kind of 
black hole solutions which are open and thus of course much more relevant in order to deal 
with the astrophysical and cosmological issues discussed at the introduction. These solutions 
are characterized (in addition to the horizon) by non-vanishing value of S(oo) which fixes the 
cosmological constant parameter by k = vS 2 (oo) /2. They have the asymptotic behavior 

B( r ) ~ kt 2 + (B 2 log(r) + B x )r + B , S(r) ~ J^- + sj^-t , (6.12) 
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where Bq, B\,B2, S\ are constants. Fig. [8]contains a graphic representation of typical solutions. 
6.3 Regular solutions 

Apart from black hole solutions discussed in the previous section, the system (15. 6p also admits 
regular solutions in r £ [0, oo] which may be viewed as gravitational solitons in this scalar- 
tensor theory. Taking advantage of the symmetry (j6.3l) . the boundary conditions can all be 
fixed to be : 

B(0) = 1 , B'(0) = , B"'(0) = , 5'(0) = , B"(oo) = 2k (6.13) 

with 5(0) as an additional input which we will allow to vary in a certain range. It turns out that 
a family of solutions exists, labelled by 5(0). The solutions behave asymptotically according 
to (16. 9h and several profiles of such solitons are presented in Fig. [9l 

Actually, this kind of solutions must exist also for the boson star system of section [3] as 
completely static self- gravitating solutions - a novelty in the CG with no analogue in standard 
GR. We have not found them at the time because we used uj in order to rescale the dimensionful 
variables. Note however that the interpretation of these static solutions is very different in both 
cases: Here it is a gravitational soliton in a scalar-tensor theory, while the same solution within 
the purely tensorial CG describes a self-gravitating scalar field, i.e. boson star, with the peculiar 
property that it is purely static. The dotted line in Fig. [13] below shows the (inertial) mass 
of such a static boson star as a function of the central value of the scalar field 5(0). The 
oscillatory curve is very similar to the one found for the usual boson stars in GR [10 } [12]. 

The striking feature about these regular solutions is that the field B deviates only a 
little from the form B(r) = 1 + Kr 2 . The numerical results indicates that the difference 
1 — (B" (oo) / B" (0)) is positive and of the order of a few percents (we checked that this is not a 
numerical artefact). As a consequence, these regular solutions are essentially characterized by 
their cosmological constant. 

If we examine the family of black holes with a fixed A and decreasing r^, it turns out that 
the maximal value b max of the parameter b increases. The numerical results then strongly 
suggest that the profile of the regular solution is approached on the interval r G]0, oo[ by the 
black holes corresponding to the second branch. The convergence cannot be extended towards 
the point at the origin because of the different condition of the metric field: B{rh) = for 
black holes, B(0) = 1 for the soliton. 
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(a) (b) 

Figure 9: Several profiles of the regular solutions of the scalar-tensor system with v = for different 
central value of the scalar field: 5(0) = 1, 5(0) = 4, 5(0) = 16. Notice the oscillations for large 5(0). 



7 Perfect Fluid Solutions in Scalar-Tensor Conformal Gravity 

Having investigated the scalar-tensor vacuum solutions and especially obtaining open solutions, 
we now proceed to couple matter sources to this system. The first is the perfect fluid with a 
polytropic equation of state. In this case we are confronted with the following set of equations: 



(rB)"" 1 

- — 1 

r B 



BS 



12 



3" / „ N 



(7.1) 



£™l*S-u* = (7.2) 



which are supplemented by the conservation law (|1.16p . 
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We have found two types of regular solutions to the equations above distinguished by 
S(oo) being either zero or non-zero. This corresponds to closed or open spacetime geometries 
respectively. For the two cases, the field B(r) satisfies the same boundary conditions as in the 
purely tensorial Conformal Gravity, namely 

B(0) = 1 , B'(0) = , B"'(0) = , B"(oo) = 2k (7.3) 

The boundary conditions on the function S(r) are different for the two solutions. Setting v = 0, 
we find solutions with the asymptotic behavior 

B(r) ~ nr 2 + B x r + B , S(r) ~ -± , P r (r) ~ ^ (7.4) 

The space-time associated with these solution is closed since the function rS(r) varies on a 
finite range. Similarly, the proper radial distance is bounded from above. The solutions of 
this type can be deformed for v > 0; however, they do not exist for large values of v. The 
function S(r) indeed develops a singularity at a finite value of r when v approaches a critical 
value v = v c . 

The solutions of the second type that we constructed are open and exist for generic non-zero 
values of the coupling constant z/; they are characterized by 5(oo) > and obey asymptotically 

B{r)~nr 2 + (B 2 Iog(r) + Bi)r + B , S(r) ~ J** + S^^- , P r (r) ~ 4 ( 7 - 5 ) 

V v r 

where Bq, B\ } B2, Si are constants. This form was checked both analytically and numerically. 
The corresponding space-time is open since rS(r) is unbounded from above. Let us point out 
two features of these solutions (i) They do not possess a regular limit for v — > 0, as seen e.g. 
from 5(oo) = \Y2k/v. (ii) Non analytical terms (log terms) appear in the asymptotic expansion 
of the fields B and S. These terms seem to be related to the fact that the field S does not go 
to zero for r — > 00. 

Fig. [TUh shows the profiles of two solutions with v = 1. The functions B(r) and P r (r) are 
quite similar to those in the "pure tensor" theory. The main difference is that point particles 
are now consistently coupled to the gravitational field through the new field S(r). The coupling 
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Figure 11: Perfect fluid solutions in the scalar-tensor theory: dependence on the parameter k. The 
other parameters are: v = 0.2, A = 1/3. 
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is described now by the combination S 2 (r)B(r) and it is obvious that besides the cosmological 
r 2 behavior, we recover the linear potential with logarithmic modifications. Further study is 
required in order to check the relation of this new kind of solutions to observational data. Fig. 
[TUb presents the dependence on the polytropic index A of the main properties of the solutions. 
These properties as a function of the cosmological constant parameter k are shown in Fig. [TTT 



8 Boson Stars in Scalar- Tensor Conformal Gravity 



Next we move to the complex scalar field, i.e. boson stars. In this case, the Lagrangian density is 
the sum of all the previous terms with a possible additional coupling between the "gravitational 
scalar field" S{x) and the other scalar <£, namely — /j,S 2 \$\ 2 /2 with fi dimensionless (real) 
parameter. The field equations will contain two second order equations for the two scalar fields 



(r 2 Bf)' 



+ 







{r 2 BS')' R 



and a fourth order equation for B(r) 
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2B\ 
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3.2 



r I 
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r 



i 



^ 2 / 2 + i(^'+v) (z 2 )'-^/ 2 



(8.3) 



The space of solutions is quite large and defined by two types of parameters: those which 
appear in the field equations namely, A and \i and parameters (integration constants) which 
specify the solutions like 5(oo), n etc. 

A systematic survey of all possible solutions is beyond the scope of this work. Here we 
present the main properties of several families of solutions in limited but typical regions of 
parameter space. Here too we find closed as well as open solutions. 

We addressed the system (j8.ip - (|8.3jl numerically. We started by fixing the different coupling 
constants according to v = [i = 0, the constant a can then be set to a = 1 by a rescaling of /. 
A rescaling of the radial variable and of the field S{x) allows one to fix u = 1. In the reduced 
system fixed this way, we further assumed A = 1. Regular solutions to the equations can then 
be constructed with the following boundary conditions : 



B(0) = 1 , B'(0) = , B"'(0) = , B"(oo) = 2k 



(8.4) 



for the metric function, and 



S(0) = So , 5'(0) = , /'(0) = , f(r oo) 



(8.5) 



for the two scalar functions. Here So is an arbitrary constant. In the numerical analysis, we 
set k = 0.1. In the limit Sq — > 0, we have S(x) = and the boson star solutions of sec. [Hare 
recovered. We have studied how the boson star solution available in CG is deformed by the 
additional scalar field and discovered a rather unexpected pattern which we now discuss. 
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S(0) 

(a) (b) 

Figure 12: Boson stars in the scalar-tensor theory corresponding to A = 1, v = fj, = 0, K = 0.1. (a) 
mass, particle number Q and the value /(0) as a function of 5(0); (b) details of the profile for S(0) — 100. 
The high value of S(0) was chosen to get noticeable oscillations - see text. 

8.1 Closed Solutions 

Increasing the parameter So gradually, we observe that the boson star gets continuously de- 
formed by the new scalar field. It turns out that the scalar field of the boson star tends 
uniformly to the null function for some critical value So = S c ; with our values of the coupling 
constants, we find S c « 6.26. Accordingly, the mass corresponding to the boson star (i.e. sup- 
ported by the field f{r)) tends to zero in this limit, along with the particle number Q. These 
features are illustrated by Fig. PT2"h . So, for < So < S c two regular solutions coexist: the 
scalar-tensor boson star (STBS) and the scalar-tensor regular solution (STR - the gravitational 
soliton) with f(r) = . The values of B"(0) and the inertial mass Mis of the scalar field S 
(only) are represented in Fig. [T3l 
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5 10 15 20 



S(0) 

Figure 13: The value S"(0) and the inertial mass of the field S for STBS (on-line red) and STR soliton 
(on-line black) as functions of S(0). Notice the gaps in the STBS (red) curves where no solutions exist. 
This can be used to distinguish between the cases in a B&W plot. 
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However this is not the end of the story. Indeed, while we continue to increase Sq, it turns 
out that STBS solutions reappear for Sq > 9.1 and that again, the two non trivial solutions 
exist for 9.1 < Sq < 13.8. After still another gap of non existence of STBS solutions they again 
reappear for So > 15 and seems then to coexist with the STR (regular) solution as suggested 
by Figs. H2i, QJ 

It is challenging to find a full analytical explanation of the features just discussed above, 
but so far we have not fully succeeded. Possibly, an explanation is to be found in the fact that, 
for large values of 5*(0) the function B(r) develops some oscillations near the origin. These 
oscillations, which appear more clearly when looking at B'(r) and B"{r) (see Fig. 112b) are due 
to a term 

S 2 

Y" + —Y + subdominant terms = , Y = B" (8.6) 
B 

contained in the equations and which leads to visible oscillations when Sq is sufficiently large. 
A detailed numerical study of the "gap-structured" phenomenon of the STBS solution shows 
that: (i) the function B"{r) is monotonically decreasing for Sq < 6.26; (ii) a local minimum 
of B" occur somewhere for Sq E [9.1, 13.8]. This strongly suggests a connection between the 
oscillations of the function B" and the pattern of STBS solutions. 

In fact, the occurrence of oscillations appears for the STR solution already, i.e. in the 
absence of boson star. This property is illustrated in Fig. [9] where the profiles of B, S, B' , B" 
are superposed for three values of S(0). 

8.2 Open Solutions 

Finally we turn to the very different kind of solutions, namely the open ones. Along with the 
case of polytropes discussed in sec. El the open solutions are characterized by S(oo) = \Jlnjv 
and the corresponding asymptotic expansion presents log-terms. The pattern of solutions of 
these non linear equations is rich and presents several bifurcations in the space of coupling 
constants. Fixing the different coupling constants, the solutions are even not unique since 
they are characterized by parameters at the boundary (or integration constant like k or S(0)), 
generating continuous families of solutions. 





(a) 



(b) 



Figure 14: Open boson star solutions in the scalar-tensor theory: (a) two profiles for v = 0.2 and v = 1 
with A = l, [i = 0, k = 0.1; (b) plots of several characteristics of the solutions as a function of the 
parameter A with v = 0.2, /i = 0, k = 0.1. 
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Figure 15: Open boson star solutions in the scalar-tensor theory: dependence on the parameter k. The 
other parameters are: A = 1, v = 0.2, fj, = 0. 

Figure [TSa shows two typical open boson star profiles for two values of v. The quite weak 
dependence on v makes it sufficient to present in Fig. 114b the dependence of the mass, charge 
and other quantities on the self-coupling A. On the other hand, the physical quantities are 
quite sensitive to the parameter k as is clearly apparent from Fig. [T5l We see that boson 
stars exist only up to a maximal value of k and that the STBS system bifurcate into a regular 
scalar-tensor solution (STR). 

9 Conclusion 

We have analyzed several types of spherically symmetric solutions in the "minimal" CG and 
in a scalar-tensor extension. The polytrope solutions in the "minimal" case have an asymp- 
totically linear gravitational potential and contain a "wrong sign" Newtonian component of 
a 1/r term only in asymptotically anti-de Sitter space which is generated by a "cosmological 
integration constant" («). These are in line with previous results showing that the generic 
exterior gravitational fields in this theory have a behavior which is very difficult to settle with 
observations. 

The linear potential may be considered an advantage in the sense that applying it in a 
galactic scale may provide an explanation for the rotation curves without invoking dark matter. 
However, since there is no a-priori reason to refrain from trying CG in a solar system scale, 
the absence of the (attractive) Newtonian potential seems to be a drawback. Of course, it is 
always possible to claim that "conformal polytropes" made of matter which satisfies Tff = is 
a rather special kind of material that may be found in galactic and intergalactic scale, but this 
line of argument effectively brings dark matter through the "back door" . 

From this point of view, a scalar field can be viewed as a more conventional matter source, 
as scalar fields are ubiquitous in theoretical physics. Since the conformal coupling to gravity 
adds to the energy-momentum tensor terms which render the energy density non positive 
definite, it might be expected that the gravitational fields of boson stars turn out to have 
different behavior. However, we found that boson stars exist only in (asymptotically) anti-de 
Sitter space, otherwise the mass and charge of the solutions do not converge. Similarly, the 
gravitational potential of these boson stars is linear with a +l/r additional contribution. 
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The motivation for the scalar-tensor extension was to search for possible different behaviors. 
Moreover, from a formal point of view the scalar-tensor CG is the simplest theory of its kind 
which has physically meaningful time-like geodesies, or in other words, couples consistently to 
point particles. 

The pattern of the solutions we have found in this case turns out to be very rich and presents 
unexpected features, especially with respect to the purely tensorial CG. First, in the vacuum 
sector we found Schwarzschild-like (or more accurately, Schwarzschild-AdS-like) solutions as 
well as closed solutions with finite radial extension. In addition, there exist regular ("soliton- 
like") solutions which have no analogue in ordinary GR. 

When matter sources are added, the resulting solutions are classified similarly for both 
perfect fluid polytropes and boson stars: There are closed solutions which although interesting 
on their own right cannot be considered as relevant in a four-dimensional astrophysical or 
cosmological context. A second type is open solutions with a gravitational potential which 
contains the "standard" (by now) linear term modified with logarithmic corrections whose 
observational relevance needs further study. 

The kind of equations we have solved (fourth-order) is unconventional but could be treated 
with a good accuracy by our numerical methods which appear in this case to be indispensable. 
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